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COMPUTERIZED SYSTEM AND METHOD FOR PREDICTING MORTALITY 
RISK USING A LYAPUNOV STABILITY CLASSIFIER 

CROSS-REFERENCE TO RELATED APPLICATIONS 
[0001] This application claims the benefit of U.S. Provisional Application No. 

60/468,765 filed May 8, 2003. 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR 

DEVELOPMENT 

[0002] Not applicable. 

TECHNICAL FIELD 

[0003] The present invention relates to a system and method for identifying which 

individual persons admitted to an intensive care unit (ICU) are likely to die prior to hospital 
discharge, and which persons are likely to recover and survive to discharge. 

BACKGROUND OF THE INVENTION 
[0004] A very difficult clinical problem facing clinicians is knowing when further 

treatment is futile and no longer appropriate in a patient who has developed severe 
complications after surgery and is being treated in an ICU. It is now possible to prolong the 
process of dying among such patients. This results in unnecessary pain and loss of dignity 
for the patient, anguish and distress for the patient's relatives and is dehumanizing for the 
clinical and nursing staff. It has also tremendous implications in the use of limited health 
care resources. 
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[0005] The complementary clinical problem facing clinicians, that of knowing 

when further aggressive treatment is appropriate and has reasonable odds of saving the 
patient's life and yielding quality of life and other benefits in addition to survival, is 
equally difficult. Prognostic criteria based on "static" analysis of group statistics are of 
little value in decisions to withhold or withdraw therapy from ICU patients too ill to 
benefit, since they do not provide adequate information on the features that distinguish 
non-survivors from survivors. 

[0006] In 1990, Smedira and co-workers (Smedira 1990) wrote that "life support is 

withheld or withdrawn from many patients especially those with critical illness, but the 
exact number is not known." The terms 'withholding of life support' and 'withdrawal of 
life support' refer to the process according to which various medical interventions either 
are withheld from patients or discontinued with the expectation that the patient will die as a 
result. 

[0007] The feeling that a patient admitted to the ICU should not be considered as a 

'terminal patient' produces difficulties associated with rationalizing the use of new and 
costly technologies in a situation where resources are scarce, such as the modern healthcare 
systems. This has led to the development of new systems for the assessment of severity for 
use in the ICU. These tools may be used to predict outcome, although their use has ethical 
and financial implications. 

[0008] Over the last 15 years, several systems for the assessment of multiorgan 

failure have been developed, mainly because multiorgan failure represents the main cause 
of mortality and morbidity among critically ill patients managed in the ICU. 
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[0009] A computer model designed to improve the process of decision-making 

concerning the continuation or escalation of aggressive ICU interventions, or of abating 
interventions that are probably futile, must possess the following properties: it must reflect 
the time-dependent, or "dynamic", pathophysiological process and be able to predict death 
with high accuracy, early in the clinical course when there is but a scanty timeseries 
comprised of very few (<10) values of predictive scores or physiologic indices of end- 
organ function. It is important to emphasize that serial scoring is not purely a reflection of 
the inherent potential of the patient's organ systems to recover function sufficient to sustain 
life; it reflects the ability of an intensive care unit to stabilize the patient or reverse the 
physiological dysfunction that is present. 

[0010] Prior art includes Chang's algorithm (Chang 1988), which used 

computerized dynamic trend analysis of daily organ failure scores (APACHE II score, 
corrected for the number and duration of organ failures), noting the rate of change in score 
relative to that of the previous day and an absolute threshold to predict death. The 
algorithm was developed by tracking the daily scores of 200 ICU patients until their death 
or discharge from the ICU. It was subsequently validated on 831 patients. Chang's 
approach had poor sensitivity (predicted only 38% of all deaths) with acceptable statistical 
specificity (no false-positive predictions). More recent trend analysis techniques such as 
SOFA ('sequential organ failure assessment', see Cook 2001, Hutchinson 2000, Pettila 
2002, Rosenberg 2002) continue to suffer from inadequate sensitivity. 

[0011] Other algorithms using serial scoring have had unacceptably high false- 

positive rates. For example, Atkinson's study (Atkinson 1994) showed a false-positive 
diagnosis rate of 4.4%. If used prospectively, this algorithm does have the potential to 
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indicate the futility of continued intensive care but at the high cost of nearly 1 in 20 
patients who would survive if intensive care were continued. 

[0012] Chang's and others' approaches were unsatisfactory in terms of excessive 

reliance on long timeseries (74% of predictions resolvable by trend analysis required data 
from seven days or more in the ICU, too long to be of significant help in the contemporary 
situation with its focus on prompt decision-making soon after admission and aggressive 
discharge-planning). In part, the failure of the prior art can be traced to inappropriate 
pooling of data from groups of patients with markedly differing mortality rates, including 
many whose probability of in-hospital death was low or moderate. For this reason and 
many others, there is a need for a system and method overcoming the deficiencies of the 
prior art. 
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SUMMARY OF THE INVENTION 

[0033] The present invention is a method and system mitigating the limitations 

enumerated above and suitable for an improved system and method for predicting mortality 
risk using a Lyapunov stability classifier. 

[0034] In one embodiment of the present invention, a method in a computing 

environment for effecting a statistical assessment of mortality-predictive patterns in 
longitudinal timeseries data from individual persons admitted to hospital-based intensive 
care is provided. The method includes the following steps: accessing mortality-predictive 
serial data received from a plurality of scores; performing spectral analysis; calculating the 
Lyapunov exponent, and if the Lyapunov exponent is negative, outputting values for the 
exponent for a plurality of times in the timeseries. 
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[0035] Additional advantages and novel features of the invention will be set forth in 

part in a description which follows, and in part will become apparent to those skilled in the 
art upon examination of the following, or may be learned by practice of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0036] The present invention is described in detail below with reference to the 

attached drawing figures, wherein: 
[0037] FIG 1 is a flow chart representative of a method for predicting mortality risk 

in accordance with an embodiment of the present inventions, namely a method using 

spectral analysis and MAPLE-based determination of Lyapunov exponents; 
[0038] FIG. 2 is a diagram of Lyapunov exponents of non-survivors in the 

development dataset, for serial APACHE score timeseries as a function of time in days 

following admission to ICU; and 
[0039] FIG. 3 is a diagram of Lyapunov exponents of survivors in the development 

dataset, for serial APACHE score timeseries as a function of time in days following 

admission to ICU. 

DETAILED DESCRIPTION OF THE INVENTION 
[0040] The present invention is a method and system mitigating the limitations 

enumerated above and suitable for an improved system and method for predicting mortality 
risk using a Lyapunov stability classifier. The method and system are preferably 
implemented in one of a variety of known general purpose or special purpose computing 
system environments or configurations. Those of ordinary skill in the art will appreciate 
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that any of a variety of components and interconnection are well known, and details 
concerning the construction need not be disclosed in connection with the present invention. 

[0041] When "multiple organ systems failure" (MOSF) is defined as severe 

physiologic abnormalities, it is found to arise in the context of a wide variety of diseases. 
The distinguishing feature of MOSF appears not to be the underlying etiology, but the 
uniform and frequently fatal outcome once it develops. This indicates that MOSF may 
represent a final common pathway to death rather than a clinical syndrome with a single 
underlying etiology. Regardless of cause, evaluating the incremental impact of new therapy 
on MOSF will be assisted by the adoption of uniform definitions as well as by the explicit 
measurement of severity of disease and calculation of individual probabilities of death. 
And, regardless of cause, the fact that the pattern of physiologic derangements are 
hallmarks of the final common pathway of MOSF enables application of a unified 
approach to timeseries analysis. 

[0042] The present invention resulted from a study of the incidence of withholding 

or withdrawing therapeutic measures in intensive care unit (ICU) patients, as well as the 
possible implications of organ failure assessment (APACHE m and SOFA) in the decision- 
making process and the ethical conflicts emerging from these measures. 

[0043] The patients (n = 372) were assigned to different groups: those surviving 

one year after ICU admission (S; n = 301), deaths at home (DH; n = 2), deaths in the 
hospital after ICU discharge (DIH; n = 13) and deaths in the ICU (DI; n = 56). The last 
group was divided into the following subgroups: two cardiovascular deaths (CVD), 20 
brain deaths (BD), 25 deaths after withholding of life support (DWH) and nine deaths after 
withdrawal of life support (DWD). 
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[0044] APACHE m, daily therapeutic intervention scoring system (TISS) and daily 

SOFA scores were good mortality predictors. The length of ICU stay in DIH (20 days) and 
in DWH (14 days) was significantly greater than in BD (5 days) or in S (7 days). The 
number of days with a maximum SOFA score was greater in DWD (5 days) than in S, BD 
or DWH (2 days). 

[0045] Daily APACHE scoring is useful when the decision to withhold or withdraw 

treatment has to be considered, especially if the established measures do not improve the 
clinical condition of the patient. Although making decisions based on the use of severity 
parameters may cause ethical problems, it may reduce the anxiety level. Additionally, it 
may help when considering the need for extraordinary measures or new investigative 
protocols for better management of resources. 

[0046] The present invention accomplishes the following objectives: 

1. To determine the incidence of withholding and withdrawing care in 
critically ill patients in a multidisciplinary adult ICU. 

2. To evaluate whether decisions made about withholding or withdrawal of 
treatment and applying a pre-established protocol could be supported by the 
information obtained from the monitoring of indicators of severity, in other 
words sepsis-related organ failure assessment SOFA, TISS, and APACHE 

m. 

3. To decrease the level of anxiety among family members and the medical 
team responsible for the patient's care, at the time when the decision is taken 
to withhold or to withdraw therapeutic management. 
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[0047] Prior to developing the method disclosed herein, deaths among patients 

admitted to the ICU were defined in relation to the type of management before death. 
[0048] The patients were placed in groups according to the final outcome: those 

surviving one year after ICU admission (S); death at home after hospital discharge (DH); 

death in the hospital once discharged from the ICU (DIH); and death in the ICU (DI). 
[0049] In the DI group, the cases were classified into the following subgroups 

according to the type of management decision taken at the time of death : 

1. Cardiovascular death (CVD): those patients who died despite 
cardiopulmonary resuscitation (CPR) for at least 30 minutes, according to 
hospital guidelines. This group represents a substantial investment in 
therapeutic resources at the time of death. 

2. Brain death (BD): those patients where an irreversible loss of cerebral and 
brainstem function were found, following the guidelines of the American 
Society of Anesthesiologists and U.S. law. Once BD was confirmed, 
mechanical ventilation was withdrawn or was continued until organ 
donation was authorized. 

3. Death after withholding of life support (DWH; patients in whom it was 
decided to limit therapy): the most common therapies withheld are CPR in 
the event of cardiac arrest and dialysis. Less frequently withheld therapeutic 
measures include mechanical ventilation, laboratory analysis, surgical 
procedures, administration of antibiotics, parenteral nutrition, blood 
transfusions, fluid therapy, antiarrhythmic drugs or vasopressors. In the 
ICU, the withholding of any therapeutic measure besides CPR and dialysis 
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is rare, because either these cases would not be considered to be candidates 
for admission into the ICU, or the withdrawal of therapeutic measures of life 
support would be considered. 
4. Death after the withdrawal of therapeutic life-support measures (DWD; 
patients in whom a gradual withdrawal of therapeutic resources of life 
support was carried out, after the previous establishment of withholding 
measures): initially, nutrition, vasopressor agents and dialysis were 
withdrawn. Finally, oxygen support was withdrawn and, if necessary, the 
respiratory frequency and the tidal volume were reduced. The patient was 
never disconnected from mechanical ventilation, nor was sedation stopped 
in any case. 

[0050] The ICU doctor responsible for the patient's care, on their own initiative or 

after considering a proposal from the patient's family or the head of the department, 
proposed the need to establish therapeutic restrictions or to withdraw treatment. This 
opinion was discussed with the other doctors of the service and the nursing staff 
responsible for the patient. If the proposal was accepted, the relevant doctors from other 
departments who had sought the patient's admission to the ICU were asked for their 
opinion and, once a consensus was reached, this was conveyed to the family by the ICU 
doctor directly responsible for the patient. 

[0051] If at any time the decision was not accepted, the required life-support 

measures were continued. The ICU doctor on call always respected the decisions reached. 
Only in rare circumstances did the on-call doctor, in agreement with the admitting doctor, 
decide, with the family, to withhold treatment before group discussion. 



CRNC.2003.008 



11 



V * 

[0052] Once the patient was discharged to a general ward, the ICU doctor informed 

the ward doctor about the decisions taken, such as not to return to the ICU in case of 
deterioration in the patient's clinical condition or not to perform CPR. 

[0053] Resource utilization was evaluated by measuring the use of mechanical 

ventilation and the daily TISS in the ICU. The daily APACHE m and the daily SOFA were 
used as indicators of severity. APACHE and SOFA scores were calculated for each organ 
system. 

[0054] All patients with an ICU stay longer than 24 hours were included. 

[0055] The groups were compared using the Kruskal-Wallis test and the Mann- 

Whitney rank-sum test, with the Tukey-B test and the Scheffe test correction for between- 
groups analysis. The qualitative variables were analyzed with the Chi-square test. 
Differences were considered significant at P < 0.05. Groups with small samples (< 5) were 
excluded from the statistical analyses. All statistical analyses were performed on a personal 
computer with SAS Version 8 and MAPLE Version 8. 

[0056] A total of 623 patients were admitted to the ICU during the study period 

from January 2000 through April 2002. Two hundred and fifty-one patients did not fulfill 
the inclusion criteria and were excluded. 

[0057] A total of 372 patients were included: 68% were male, mean age (± SD) was 

59 ± 17 years (range 14-92), mean average stay in ICU was 8 ± 10 days (range 1-81) and 
the 1-year mortality rate was 59% (71 cases). The mortality rate was also closely related to 
the need for mechanical ventilation. Of 203 cases (55%) who required mechanical 
ventilation, 106 died (52%). 
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[0058] In patients who died, APACHE EI (BD: 91 ±31, range 34-129; DWH: 84 ± 

3, range 31-153; DWD: 81 ± 25, range 41-113; DM: 91 ± 31, range 44-137) and SOFA 
score on the first day (BD: 9 ± 3, range 2-16; DWH: 9 ± 4, range 1-17; DWD: 9 ± 4, range 
1-16; DIH: 10 ± 4, range 3-16) were significantly higher than in the S group (APACHE 
DI: 38 ± 24, range 2-117; SOFA score on the first day: 4 ± 4, range 0-14; P < 0.0001). The 
CVD group (APACHE D3: 41 ± 19, range 27-54; SOFA score on the first day: 2 ± 3, range 
0-4) had similar results to the S group. 

[0059] The TISS score on the first day is significantly higher in the BD (34 ± 10, 

range 15-53), DWH (34 ± 11, range 16-57) and DIH groups (34 ± 5, range 25-41) than in 
the S group (22 ± 13, range 4-54; P < 0.0001). There were no significant differences 
between the S group and the DWD (24 ± 13, range 8-46) and CVD groups (15 ± 1, range 
14-16). 

[0060] The mean SOFA of the S group (3 ± 3, range 0-12) and the CVD group (4 ± 

5, range 0-7) was significantly lower than in the other groups (BD: 11 ± 3, range 5-19; 
DWH: 12 ± 4, range 4-18; DWD: 11 ± 4, range 7-19; DIH: 8 ± 3, range 3-11; P < 
0.0001). 

[0061] Similarly the mean TISS of the S group (21 ± 1 1, range 4-47) and the CVD 

group (21 ± 14, range 1 1-31) was significantly lower than in the other groups (BD: 35 ± 9, 
range 22-60; DWH: 37 ± 11, range 15-59; DWD: 33 ± 7, range 19-44; DIH: 31 ± 5, range 
25-41; P< 0.0001). 

[0062] ICU stay (days) was significantly shorter in the S (7 ± 8, range 1-58), CVD 

(4 ± 3, range 2-6) and BD groups (5 ± 6, range 1-20), compared with the DWH (14 ± 16, 
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range 1-81) and DM groups (20 ± 14, range 5^44; P < 0.0001). There were no significant 
statistical differences between the S and DWD groups (12 ± 10, range 3-29). 

[0063] BD (30%) and DWD (15%) were differentiated from one another, given that 

the medical care and ethical implications are very different in these two groups. The DWH 
group represented 35% of cases. While brain death caused only occasional diagnostic 
difficulties, the greatest ethical conflicts among the medical team responsible for the 
patient were caused when deciding in which patients and at what moment should 
therapeutic restrictions or withdrawal of life-supporting measures be suggested and 
established. In our study, these ethical discussions arose in 48% of patients who died. 

[0064] At this point, it would be easy to generalize decision making according to 

these markers, as they seem to be good indicators of bad prognosis. However, decisions 
should be individualized for each patient, as the concepts of quality of life and dignity in 
death, new technologies, greater demand for health or the maintenance of sustainable or 
equitable medicine greatly influence the decision-making process. 

[0065] With a follow-up based on objective and quantifiable data, the doctor is 

better informed and more accurate decisions are achieved. This is fundamental when the 
changes in the role of the doctor in the decision-making process are considered. The 
principles of autonomy and self-determination, and the new concepts of quality of life, lead 
to consensual medical treatment and to adequate informed written consents. 

[0066] However, an APACHE score alone, or even trend analysis of serial 

APACHE scores, does not provide sufficient complementary information about the 
possible need to withdraw treatment, in particular when referring to the withdrawal of life 
support. Beginning a treatment is easy, but stopping it is much more difficult. The only 
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analyzed data which may be of interest in this sense are the number of shifts with a high or 
increasing score. This has high sensitivity as a sign of diagnostic or treatment failure if 
accompanied by a Lyapunov exponent indicating instability. 

[0067] In this sense, although it may be very difficult to accept a feeling of failure, 

the data afforded by Lyapunov assessment of serial scores may be very useful. The 
underlying study suggested that a patient suffering from the failure of six organs together 
with unfavorable APACHE scores inexorably evolved to death. These data were reinforced 
when the increase in the patient's APACHE score was accompanied by an equally 
significant increase in the TISS. The worsening of the clinical situation was maintained 
despite the additional therapeutic resources used. 

[0068] Undoubtedly, the withdrawal of life support reduces both ICU stay and the 

investment of therapeutic resources (lower accumulated TISS), which has corresponding 
economic repercussions. However, this reduction was not statistically significant among 
the patients in the DWH and DIH groups. 

[0069] Analysis of the data showed that patients in the DWD group were in fact 

subject to therapeutic restrictions from the moment of admission. This can be explained by 
the fact that, with similar levels of severity (APACHE EQ and SOFA score on the first day 
of ICU stay), the TISS score on the first day was lower than in the other groups of deaths. 
Despite this, the mean TISS was very similar. It was surprising to note that, during ICU 
stay, an increase of resources used was observed for the DWD group which, in the end, 
reached the mean for the other groups of deaths. 

[0070] Another noteworthy finding was the length of time that elapsed from the 

beginning of therapeutic restrictions until the decision to withhold life support. It was 
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shown that patients in whom life-support measures were withheld had previously passed an 
average of four days with constantly high score levels. This may be seen as a worryingly 
long period when costs and the optimal use of scarce resources are considered. It should 
also not be forgotten that, during this period, the patient is dying. It is not easy to define the 
'best moment' to withdraw life-support measures without overextending an 'agonizing 
state' that leads to discomfort for the patient or the family. The withdrawal of all life- 
support measures except for mechanical ventilation may mean that the patient is 
maintained in a terminal situation for some time. If this situation is maintained for a few 
minutes or even hours, it may facilitate contact between the patient and the family during 
the last moments, but if prolonged for more than a night or for days, it generates significant 
anxiety among both the family and the team caring for the patient (Christakis 1993). The 
present invention enables this difficult decision-making to be approached with objectivity 
and greater confidence than has heretofore been practical with prior art scoring systems and 
trend analysis. 

[0071] Lyapunov stability analysis is well known in applications involving 

nonlinear control systems, including military defense systems and nuclear reactor control 
systems (De Queiroz 2000, Freeman 1996, Lakshmikantham 1991, Pai 1981, Wihstutz 
1986), but it has not previously achieved practical application in analysis of medical 
variables or clinical decision-support. Exploratory consideration has appeared, as related 
to Electroneurodiagnostic (EEG) analysis (Elger 2000) and electrocardiography (EKG) 
analysis (Goldberger 1996). The Lyapunov exponent is a quantitative measure of 
separation the trajectories that diverge widely from their initial positions. The magnitude of 
the Lyapunov exponent is related to how chaotic a system is. The larger the exponent, the 
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more chaotic the system. For periodic signals, the Lyapunov exponent is zero. A random 
but stable signal will also have an exponent very close to zero. 

[0072] Spectral analysis alone cannot distinguish a chaotic process either, but some 

investigators have suggested that a particular spectral pattern (one in which the power 
density is inversely related to frequency) is highly suggestive of a nonlinear or chaotic 
process (Goldberger & West 1987, Goldberger 1996). However, the diagnostic value of 
this 1/f pattern has also been questioned (Pool 1989). 

[0073] Referring to FIG 1, a flow diagram is provided which illustrates an 

embodiment of a system and method for generating the Lyapunov exponent classifier and 
verifying and validating whether such a detector achieves statistical sensitivity and 
specificity in the intended mortality range of deployment, sufficient for satisfactory 
performance in the use for classifying patients according to in-hospital mortality outcome. 

[0074] In a preferred embodiment, a system and method of the invention includes 

the following steps: 

[0075] Let L(xi, X2, .., x n ) be a scalar function of n components of x. L(x) is 

positive-definite in a neighborhood N of the origin if L(x) > 0 for all x £ 0 in N and L(0) = 
0. 

[0076] Let x*(t) = 0, t > to be the zero solution of the homogeneous system x0 = Ax 

where x(0) = x 0 = 0. Then x*(t) is globally stable for t > to if there exists L(x) with the 
following properties in some neighborhood N of 0: 

(i) L(x) and its partial derivatives are continuous. 

(ii) L(x) is positive-definite, or L(x) > 0. 

(iii) dL(x)/dt is negative-definite, or dL(x)/dt < 0. 
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[0077] By (ii) the quadratic form L(x) exhibits an ellipsoid curve. By (iii), the 

ellipsoid curve shrinks to zero. Choose e > 0 such that Ne c N above. Any half-path 
starting in Ne remains in it because L(x) is a quadratic form (by (ii)) which exhibits an 
ellipsoid curve that is continuous as well as its partial derivatives (by (i)). The same holds 
for every sufficiently small e > 0 and hence for every sufficiently small neighborhood of 
the origin. The zero solution is therefore globally stable. This is the case for APACHE HI 
score timeseries of 99% of patients who are survivors to discharge from the hospital. 
Lyapunov exponents H for patients who did not survive were uniformly negative at some 
point during the ICU stay, and serve as a leading indicator of downward physiologic spiral 
toward death. 

[0078] In other words, the system (dx/dt) = Ax is globally stable if and only if for 

some positive-definite matrix W, the equation: 

A* H + HA = -W 

has a positive-definite matrix H. 
[0079] If for some positive-definite matrix W, the equation A* H + HA = -W has a 

positive-definite matrix H, let us show that (dx/dt) = Ax is globally stable. Since H is 
positive-definite, then L(x) = x x Hx is positive-definite (where x l is now the transpose of x 
and not the time derivative), i.e. L(x) > 0. Also, L(x) positive-definite implies that V(x) and 
its partial derivatives are continuous. Differentiating L(x), then: 

dL(x)/dt = (dx l /dt)Hx + x' H(dx/dt) 

or, as dx/dt = Ax: 

dL(x)/dt = (Ax) 1 Hx + x l HAx 
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= x l A 1 Hx + x l HAx 

= x l (A 1 H + HA) x 
thus, as A 1 H + HA = -W: 
dL(x)/dt = x l (-W)x 

W determined to be positive-definite implies that -W is negative-definite, thus: 

dL(x)/dt = x l (-W)x < 0 

[0080] Finally, it is notable that (i) L(x) and its partial derivatives are continuous; 

(ii) V(x) is positive-definite; (iii) dL(x)/dt is negative-definite. As a result, dx/dt is globally 
stable according to our previous theorem. Conversely, if dx/dt = Ax is stable, then for some 
positive-definite matrix W, the equation A 1 H + HA = -W has a positive-definite matrix H. 
dx/dt = Ax stable implies all the eigenvalues of A are negative, i.e. X < 0 for any 
eigenvalue X of A. Now, as X x = Ax, then (Ax) 1 = (k x) 1 , which implies x ! A 1 = X x l . Thus, 
premultiplying A 1 H + AH by x* and postmultiplying it by x, the following is obtained: 

x l (A 1 H + HA)x = x l (-W)x 

or: 

x l A 1 Hx + x l HAx = x l (-W)x 
or substituting in X x ! and X x: 
X x' Hx + x l HX x = x l (-W)x 
or simply: 
2X x l Hx = x l (-W)x 
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[0081] As -W is negative-definite, then x l (-W)x < 0, thus 2\ x l Hx < 0. As X < 0 by 

the assumption of stability, then it must be that x l Hx > 0, or H is a positive-definite matrix. 

[0082] Thus, we have proven that a real n x n matrix A is a stable matrix if and 

only if there exists a symmetric positive-definite matrix H such that A 1 H + HA is negative- 
definite. In practice, a choice of W = I may be made and H can be solved and solve for H in 
the equation A 1 H + HA = -I. The solution has the form H = a (A 1 )' l A~ l + P I where a and 
p are constants. Thus, choosing a Lyapunov function, L(x) = x l Hx, this solution is used to 
determine H. 

[0083] A preferred embodiment of the present invention utilizes a source code 

implemented in the commercially available mathematics software package MAPLE. 
Practical implementation as part of an integrated ICU information system entails passing 
parameters and timeseries information to the Lyapunov program, which is executed by a 
server running MAPLE in client-server or batch mode. 

[0084] A preferred embodiment of the present invention is attached as source code 

to provide an exemplary means of implementing the present invention. In the preferred 
embodiment, the second-order polynomial function f(x) = r*x*(l-x) is utilized to represent 
the system whose stability is characterized by the present invention. Depending on the 
system, the system may be characterized by a function of different order or form. If the 
structure of a particular system is not known, the structure may be developed by Taylor 
series regression, spectral analysis or timeseries analysis techniques or other methods of 
modeling known to those of skill in the art. 

[0085] Although the invention has been described with reference to the preferred 

embodiment illustrated in the attached drawing figures, it is noted that substitutions may be 
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made and equivalents employed herein without departing from the scope of the invention 
as recited in the claims. For example, additional steps may be added and steps omitted 
without departing from the scope of this invention. Also, it will be understood by those 
skilled in the art that other scores such as the Simplified Acute Physiology Score (SAPS), 
SAPS II, or MPS may be used without departing from the scope of the invention. 
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